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Using the formal languages Schoonschip and Form, we have developed general codes that 
are able to carry out all the algebraic manipulations needed to perform analytic lattice 
calculations, starting from the elementary building blocks (propagators and vertices) of 
each Feynman diagram. The main difficulty resides in the fact that, although there are 
many built in instructions to deal with Dirac gamma-matrices, Schoonschip and Form 
have been conceived having in mind a continuum theory, which is invariant with respect 
to the Lorentz group. On the lattice, on the contrary, a field theory is only invariant 
with respect to the hypercubic group, contained in the (euclidean) Lorentz group and not 
every pair of equal indices should be summed over. Being impossible to directly use the 
'gammatrics' of Schoonschip and Form as they are, special routines have been developed 
to correctly treat gamma matrices on the lattice, while using as much as possible of the 
built in Schoonschip and Form commands. 

We have used our codes to compute, in 1-loop perturbation theory in lattice QCD, 
the renormalization constants and mixing coefficients of the operators that enter in the 
determination of the first two moments of deep inelastic scattering structure functions. 

1. Introduction 

In this contribution we present a short report on the use of some algebraic manip- 
ulation computer languages, namely Schoonschip and Form, to carry out the large 
and complex manipulations which occur in certain analytic perturbative lattice cal- 
culations. We have employed our codes to compute in lattice QCD matrix elements 
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of operators related to the first two moments of the structure functions measured 
in Deep Inelastic Scattering (DIS), both using the standard Wilson action and the 
improved nearest-neighbor action. Up to now we have completed the calculations 
of the 1-loop renormalization constants and mixing coefficients of the fermion and 
gluon rank two operators (O 9 and O 9 ) that are related to the first moment of 
DIS structure functions, and of the rank three operator (0^„ r ) related to the second 
moment. 

Consequently to the use of improvement in lattice QCD, there is indeed a large 
number of diagrams to be computed and, in most of them, even in the simplest 
case of the rank two operator, algebraic manipulations are very complicated, each 
Feynman diagram giving rise to a huge number of terms. For this reason, we 
have found unavoidable to check independently, by means of suitable algebraic 
manipulation programs, all the computations made by hand. For OX„ T we almost 
completely have to rely on computer results. 

The first code we developed was written using the Schoonschip language, and 
was used in the case of O 9 and O 9 .0 Afterwards, we upgraded this code to make 
it able to perform the more complex calculations concerning 0*„ T . In this case 
it is an impossible task to calculate all diagrams by hand; the agreement for the 
few of them that we could really compute by hand was nevertheless completefl 
Further checks were later offered by an independent code developed by Beccarini 
using Formjl and by another Form code constructed by us, that we have further 
upgraded to meet the requirements of the higher complexity of the operators < j 1VT 
that we are presently studying.cl 

2. Moments of Structure Functions and Lattice QCD 

In QCD the lowest twist operators appearing in the light-cone expansion of the 
product of two hadronic weak currents, relevant to DIS, are 

°ll% N = ^v^7[ M1 IW" D MN ] (1 ± 75) y<A 
gS = Y Tr 

P 

where the A^'s are flavor matrices Bl^ttlA The operators (|l|) are gauge invariant and 
all have twist twojlJ they are symmetrized with respect to all Lorentz indices. S 
and NS superscripts refer to Singlet and Non Singlet flavor structures. 

The operators of rank two and three in (Q) are related to the first two moments of 
the DIS structure functions via the Wilson operator expansion of the product of two 
electromagnetic currents (we do not consider in this talk the 75 contributions 
The knowledge of the hadronic matrix elements of these operators is necessary for 
the theoretical evaluation of the moments of the x-distributions of quarks and gluons 
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inside the hadrons (x is the momentum fraction carried by the struck parton inside 
the hadron).E3 

Our work consists in the determination of the 1-loop renormalization constants 
and mixing coefficients of the lattice fermion and gluon operators of rank twop 
and the renormalization constant of the fermion operator of rank three that cor- 
respond to the continuum expressions given in (]!]).□ The calculations were done 
using the nearest neighbor improved lattice QCD action and include the results for 
the renormaliz ation constants and mixing coefficient pertinent to the standard Wil- 
son actionJiS'U'Ej'Ej'O'Ea All these constants are needed to renormalize the lattice 
operators and be able to extract physical hadronic matrix elements from numbers 
obtained in Monte Carlo QCD simulations. 

Lattice QCD represents today the only viable way of evaluating from first prin- 
ciples the hadronic matrix elements needed for the computation of the moments of 
the structure functions. The nearest neighbor improved QCD action (also known 
as the "clover-leaf" action), is obtained by adding to the standard Wilson action 
(for one flavor / on a euclidean lattice) 
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Tr [Un^Un+n^Un+^pU+v] — N c 



the nearest-neighbor interaction term 



Here F n is not the naive lattice "plaquette" 
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with U n ^ v = Un^Un+^vU^+v JJ^v-i but rather the average of the four plaquettes 



lying in the plane /iz/, stemming from the point n: 
F, 



(5) 



flf— ± 



liiy— ± 



The use of this action has been proved to remove from on-shell hadronic matrix 
elements all terms that in the continuum limit are effectively of order aJ13'E3iljO 
a being the lattice spacing, provided that in the calculation of a fermionic Green 
function, each fermion field undergoes the rotation 
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With these requirements, the difference between continuum and lattice is lowered 
from 
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p' ) is the physical matrix element we want to extract from Monte 

phys. 

Carlo data and d is its physical dimension) to 
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Using this recipe, the systematic error related to the latti£e_discretization drops in 
many cases from 20 - 30 percent down to 5 - 10 percent£2li23 

The numbers we want to compute are the renormalization constants that connect 
the bare lattice operators on the lattice, 0(a), to finite operators, O(fi), renormal- 
ized at a scale fi: 

&{^)^Z lk {^a)O k {a). (9) 

The constants Zi^ are fixed in perturbation theory by the same renormalization 
conditions used in the continuum. In the flavor Singlet case there is a mixing be- 
tween quark and gluon operators of rank two that have the same conserved quantum 
numbers. However, with the choice /u, ^ v and the definition (g) of the gauge field 
strength we avoid mixing with operators of lower dimensions, and therefore the need 
for subtractions of power divergences. We thus write: 



O q 



Z qq O" 



z qg o 9 



O 9 = Z gq O q + Z gg 9 , (10) 
where the Z's are determined by imposing the renormalization conditions: 



<q\Oi(ti)\q> 

< g,a\d q (n)\g,(j > 

<q\d*(ri\q> 

< gi a\d 9 (n)\g,a > 



= <q\O q {a)\q>\ t ; 2 % 2 

= <g,<j\O c >(a)\g,a>\ t ; 2 % l2 =0 

= <q\0°(a)\q> |£ e => =0 

= <g,<j\0 9 (a)\g,a> 1^2, 



(11) 



where \q > and \g, a > are respectively quark and gluon states. 

Similarly in the case of 0^„ T there is no mixing with lower dimensional operators, 
if one chooses fi v =/= t. In fact in the "quenched" approximation (no internal 
quark loops) and Ot VT are simply multiplicatively renormalizablc (no mixing 
with other operators). 



3. What our Programs do 

We have developed general codes able to automatically carry out all the necessary 
algebraic manipulations, starting from the elementary building blocks of the calcula- 
tion represented by the expressions of propagators, vertices and Fourier transforms 
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of the operators. The final output of our codes can be cast either in the form 
of an analytic expression (to be possibly compared with the calculations made by 
hand) or in a format ready to be introduced in a suitable Fortran code for the final 
numerical loop integration. 

The main difficulty in using Schoonschip and Form in lattice calculations re- 
sides in the fact that, although these languages have many built in instructions 
to deal with gamma matrices, they have been conceived having in mind a contin- 
uum theory, which is invariant with respect to the (euclidean) Lorentz group 0(4). 
On the lattice, on the contrary, where the theory is only invariant with respect to 
the hypercubic group ff(4),EJE-3 we cannot use the usual summation conventions 
and the standard rules to rearrange gamma matrices. To make an example, one 
simple and frequently occurring term, like ^2 x "/\p\smk\, is not properly handled 
by these languages. Actually, this term is wrongly reduced by Schoonschip and 
Form to ^smk\, because the first two equal indices that are encountered are by 
default assumed to be contracted. It turns out that most of the terms that arise on 
the lattice are improperly handled by a straight use of the Schoonschip and Form 
commands. 

Being impossible to directly use the "gammatrics" of Schoonschip and Form as 
they are, it has been necessary to develop special routines to correctly treat gamma 
matrices on the lattice, while using as much as possible of the built in Schoonschip 
and Form commands. From our efforts at least three sets of routines have grown 
up, and we have checked for each of them that the rather complicated Dirac algebra 
(we have products of up to seven gamma matrices) is carried out correctly. 

The CPU times needed to run our programs has greatly shrunk with the evolu- 
tion of the various codes developed in the last three years. The first Schoonschip 
code we used in the case of the rank two operator was running on a Sun 3 work- 
station, with CPU times that varied, depending on the complexity of the diagrams, 
from 20 seconds up to 5 minutes for the most complicated cases.0 This was quite 
reasonable for our purposes, but when turning to the rank three operator the time 
required byjjhe appropriately modified Schoonschip code was one order of magni- 
tude higher J3 More recently we have developed new codes written using Form that 
run on VAX- VMS machines. These programs are considerably faster, and the typ- 
ical CPU times are back in the minute range. There has been a final reduction of 
these values when we have turned to an HP-UX 9000/735 machinefl 

Another big problem we have encountered was the limitation on the working 
memory allowed by Schoonschip and Form on the different machines we have used. 
This is of primary relevance in the first stages of the computation, when all vertices, 
propagators etc. are expanded, up to the first and second order in the lattice spacing 
a, for the rank two and three operators respectively. In our case, when dealing with 
the rank three case, we have several products of up to ten trigonometric functions, 
each one to be expanded to second order in a. Their product thus contains tenths 
of thousands of monomials. A large part of them does not contribute to the final 
expression, to the order in a we are interested in, and has to be killed at the earliest 
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possible stage of the calculation. 

As a last remark, we want to say that, besides being more portable on different 
machines, Form appears to be superior to Schoonschip, at least for our kind of 
calculations (not involving 75 and limited to four dimensions). Form seems indeed 
to be faster and to allow an easier localization of the errors. We have also found 
very useful the larger set of wildcards that Form owns compared to Schoonschip. 
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